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^^ , The phenomenological model that all quark and lepton mass matrices have 
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the same zero texture, namely their (1,1), (1,3) and (3,1) components are 

Qh| zeros, is discussed in the context of SO(10) Grand Unified Theories (GUTs). 

P-t. 

The mass matrices of type I for quarks are consistent with the experimental 
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data in the quark sector. For the lepton sector, consistent fitting to the data 

of neutrino oscillation experiments force us to use the mass matrix for the 

charged leptons which is slightly deviated from type I. Given quark masses 

and charged lepton masses, the model includes 19 free parameters, whereas the 

SO(10) GUTs gives 16 constrained equations. Changing the remaining three 

parameters freely, we can fit all the entries of the CKM quark mixing matrix 

and the MNS lepton mixing matrix, and three neutrino masses consistently 

with the present experimental data. 
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I. INTRODUCTION 

The Downward and upward discrepancy in the atmospheric neutrino deficit in Super 
Kamiokande [[l| together with other neutrino oscillation experiments such as solar neutrino 
0, reactor || and accelerator j| experiments drives us to the definite conclusion that neu- 
trinos have masses. These experiments enable us to get a glimpse of high energy physics 
beyond the Standard Model. In these situations our strategy is as follows. First, we search 
for the most suitable phenomenological quark and lepton mass matrices which satisfies mis- 
cellaneous experiments in the hadron and electro- weak interactions. Next, in order to search 
for its uniqueness and for its physical implications, such mass matrices are incorporated into 
the Grand Unified Theories (GUTs). Of course phenomenological mass matrices and GUTs 
are closely correlated and the real model building is performed going back and forth between 
these two approaches. Indeed, we consider the seesaw mechanism in neutrino mass matrix 
||, which supports minimally SO(10) GUTs. Conversely SO(10) GUTs prefer the mass 
matrices reflecting some similarity in the quark and lepton sectors. In the seminal work of 
phenomenological quark mass matrix models 0, Fritzsch proposed a symmetric or hermi- 
tian matrices later called a six texture zero model which has vanishing (1,1), (1,3), (3,1) 
and (2,2) components in both the mass matrices, M u for up-type quarks (u, c, t) and M^ for 
down-type quarks (d, s, b). Here n texture zero means that two types of quark mass matrices 
have totally n zeros in the upper half of hermitian mass matrices, in this case (1,1), (1,3) 
and (2,2) in each mass matrix. However, this model failed to predict a large top quark mass. 
Symmetric or hermitian six and five texture zero models were systematically discussed by 
Ramond et.al. 0. They found that the hermitian M u and M^ compatible with experiments 
can have at most five texture zero. Before the work of Ramond et.al. nonsymmetric or non- 
hermitian six texture zero quark mass matrices model (nearest-neighber interaction (NNI) 
model) was proposed by Branco-Lavoura-Mota @, and Takasugi showed that, by rebasing 
and rephasing of weak bases, always one of M u and Md can have the symmetric Fritzsch 
form and the other does NNI form ||. Demanding to deal with the quark and lepton mass 



matrices on the same footing, we have proposed a four texture zero model ||10|| , in which 
all the quark and lepton mass matrices, M u , M^, M e and M v are hermitian and have the 
same textures. Here M v and M e are mass matrices of neutrinos (z/ e , z/^, v T ) and charged 
leptons (e,yU,r), respectively. Namely their (1,1), (1,3) and (3,1) components are zeros and 
the others are nonzero valued. This model was also discussed by Du and Xing |TTJ, by 



Fritzsch and Xing [|T^] , by Kang and Kang []T3[ , by Kang, Kang, Kim, and Kim [Q , and by 



Chkareuli and Froggatt ||15| , mainly in the quark sector. This model is compatible with the 



large top quark mass, the small quark mixing angles, and the large v^-v T neutrino mixing 
angles via the seesaw mechanism. In this article, we discuss the above four texture zero 
model embedding in the SO(10) GUTs. The SO(10) GUTs impose some further constraints 
on the mass matrices. Using those constraints we predict all the entries of the lepton mixing 
matrix and neutrino masses, which are consistent with the experimental data, in terms of 
three free parameters left in the model. 

This article is organized as follows. In section 2 we review four texture zero model. In 
section 3 we present a mass matrix model motivated by SO(10) GUTs . This model is 
combined with the four texture zero ansatzae in section 4. Section 5 is devoted to summary. 

II. FOUR TEXTURE ZERO QUARK-LEPTON MASS MATRICES 
Phenomenological quark mass matrices have been discussed from various points of view 



|Tq| . In this section we review our quark and lepton mass matrix model |TI| . The mass term 
in the Lagrangian is given by 



->M 



QR,iMuijqtj ~ qiiMdijqij - l R ,iM eij l L>j - v R>i M Dij v Lij 



^L,i) c M Ll3 v Lj - - 2 {v' R ^M Rl3 v' R ^ + H.c. (1) 



with 



Ql,R 



( \ 

u 



\ f J 



Ql,r 



L,R 



'^ 



\»J 



l L,R 



L,R 



( \ 

e 



/i 



v r / 



vl 



L,R 



( \ 



v„ 



K Vr, 



f.A 



V„ 



VV 



R 



(2) 



where M u , Ma, M e , Md, Ml, and Mr are the mass matrices for up quarks, down quarks, 
charged leptons, Dirac neutrinos, left-handed Majorana neutrinos, and right-handed Ma- 
jorana neutrinos, respectively. The mass matrix of light Majorana neutrinos M v is given 
by 



M V = M L - MlM^M D , 



(3) 



which is constructed via the seesaw mechanism || from the block-diagonalization of neutrino 
mass matrix, 



M L Ml^ 



(4) 
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We put a ansatz that the mass matrices M u ,Md,M e and M v are hermitian and have 
the same textures. Our model is different from the Fritzsch model in the sense that (2,2) 
components are not zeros and that our model deals with the quark and lepton mass matrices 
on the same footing. The mass matrices Mjj, Ml, and Mr are, furthermore, assumed to 



have the same zero texture as M v . This ansatz restricts the texture forms |T(J and we choose 
the following our texture because it is most closely related with the NNI form ||. 
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The nonvanishing (2,2) component distinguishes our form from NNI's. Thus the quark and 
lepton mass matrices are described as follows. 
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where P d = diag(e iai , e ia \ e ia ' A ), «„• = a* - o^-, and P e = diag(e i ^ 1 , e^ 2 , e i/33 ), /^ = A - /%. 

Let us discuss the relations between the following texture's components of mass matrix 
M: 



( A ^ 
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and its eigenmass TOj. They satisfy 



B + D = mi + m 2 + m 3 , 
BD — C 2 — A 2 = m\Vfi2 + m 2 m 3 + m 3 m\, 



DA 2 



-17li17l2l7l 3 . 
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Therefore, the mass matrix is classified into two types by choosing B and D as follows: 

[type I] B — 1712, D = m 3 + mi 

[type II] B = mi, D = m 3 + m 2 (9) 



In the previous paper |10j we showed that type I is compatible with the experimental data 
both for the quark and lepton mass matrices. So, we concentrate ourselves on the type I 
case. 

In the type I case (B = rri2, D = 7773 + mi), the other A and C take the following value 
from Eq.(H) 
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Transforming m\ into —m\ by rephasing, the mass matrix M becomes 



(10) 



M 







mi momo 







rn, l ?Ti 2 Tn 3 
mo — m-± 



m 2 



m 1 m 3 (m 3 -m 2 -m 1 ) 
Trio — m-^ 



rn l m 3 ( m 3 — rri 2~ Tn l ) 
ma — m-^ 



777 3 — 777! 



v /777i7772 

% /777i777 2 777 2 ^77717773 





y/mirnz 7773 — 7771 

(for 7773 ^> 7772 ^> TTlx). 



(11) 



The orthogonal matrix O which diagonalize M in Eq.(p]) as 
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The mass matrices for quarks and charged leptons, Md, M u , and M e are considered to be of 
this type I and are given by 
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Those M d , M u , and M e are, respectively, diagonalized by matrices P d O d , O u , and P e O e . 
Here the orthogonal matrices O d , O u and O e which diagonalize P d M d P d , M u , and P\M e P e 
are obtained from Eq. (|13|) by replacing mi, m 2 , m 3 by m d , m s , m b , by m u , m c , m t , and 
by m e , m M , m r , respectively. In this case, the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix V can be written as 
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where the P^ 1 factor is included to put V in the form with diagonal elements real to 
a good approximation. Furthermore, the P~ x and P q = diag(e ?<?il , e^ 2 , e 1 ^ 3 ) with 0i — 
(f) 2 = aig(Pj l O^P d O d ) 12 and 0i — 3 = arg(P d -1 0^P d d ) 2 3 come from the choice of phase 
convention as Eq. (|15"|). The explicit forms of the components of V are obtained [I0|| as 
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The lepton mixing matrix U [hereafter we call it the Maki-Nakagawa-Sakata (MNS) 



mixing matrix [17]], is given by 



( U n U 12 U 13 ^ 



U = P\0 T e P e O v 
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where the PJ factor is included to put U in the form with diagonal elements real to a good 
approximation. Here the 0„ is the orthogonal matrix which diagonalizes the light Majorana 
neutrino mass matrices M v given by Eq.(^). 

III. MASS MATRICES IN THE CONTEXT OF SO(10) GUTS 

Even if we succeeded in constructing the quark mass matrices M u and Md consistent with 
experiments, we have infinitely many mass matrices equivalent to the M u and Md which are 
defined as 

M' U = F^M U G U M' d = F^M d G d , (18) 

with arbitrary unitary matrices F, G u , and Gd in the standard SUl(2) x C/y(l) model, and 
with G u = Gd in the SUl(2) x SUr(2) x Uy(l) model. The fact that quark and lepton mass 
matrices have the same form strongly suggests that the quarks and leptons belong to the 
same multiplets. So in this section we try to incorporate our mass matrix in the context of 
SO(10) GUTs. We consider two SO(10) symmetry breaking patterns. 

(i) 50(10) -► SU{4) x SU L (2) x SU R {2) -► SU C (3) x SU L {2) x SU R {2) x £7(1) -» G„ 
(ii) 50(10) -^ SU(5) -> G„ (19) 

where G s = SU C (3) x SU L {2) x U{1). 

A. The case of 50(10) breaking down to SU(4) x SU L (2) x SU R (2) 



Here we consider the charge-conjugation-conserving (CCC) version |TJ| |T9| [^] of the 
SO (10) model in which Left-Right discrete (not manifest) symmetry is imposed. 



In the SO(10) model [21] |22| ||23| , the left- (right-) handed fermions i^L(R)i m a given i-th 
generation are assigned to a single irreducible 16. Since 16 x 16 = 10g + 120a + 126s, the 
fermion masses are generated when the Higgs fields of 10, and 120, and 126 dimensional 
SO(10) representation (denoted by 0io, 0i2o, and (fine, respectively) develop nonvanishing 
expectation values. Their decomposition under S77(4) x SUl{2) x SU r (2) are given by 



10 = (6, 1,1) + (1,2, 2), 



120 = (15, 2, 2) + (6, 3, 1) + (6, 1, 3) + (1, 2, 2) + (20, 1, 1) 
126 = (10, 3, 1) + (TO, 1, 3) + (15, 2, 2) + (6, 1, 1). 



(20) 



On the other hand, the fermion field of 16-dimensinal SO(10) representation is decomposed 



as 



16 = (4, 2,1) + (4, 1,2). 



(21) 



With respect to SU{4) x SUl(2) x SU r (2), the left - and right- handed quarks and leptons 
of a given i-th generation are assigned as 



/ \ 

U r U y Ub V e 

\ d r dy dj, e , 



L(R)1, 



(22) 



L(R) 



Fl(r)2 and F^ R ^ are likewise defined for the 2nd and 3rd generations. Note that their 
transformation properties are Fn = (4, 2, 1) and F Ri = (4, 1, 2) and that (Fu + F Ri ) yields 
the 16 of SO(10). Since (4,2,1) x (4,1,2) = (15,2,2) + (1,2,2), the Dirac masses for 
quarks and leptons are generated when neutral components in a (1,2,2) multiplet in 0io, 
(1,2,2) and (15,2,2) in 12O , and (15,2,2) in 126 of SU{A) x SU L {2) x SU R (2) C 5O(10) 
develop nonvanishing expectation values. On the other hand, the (10, 3, 1) and (10, 1, 3) in 
0126 are responsible for the left- and the right- handed Majorana neutrino masses through 
the Higgs-lepton-lepton interactions (10, 3, 1)(4, 2, 1)(4, 2, 1) and (10, 1, 3) (4, 1, 2) (4, 1, 2), 
respectively. Here the (10, 3, 1) is the SUl{2) triplet Higgs (denoted by 0(10, 3, 1)) and the 
(TO, 3, 1) is the SU R {2) triplet Higgs (0(10, 1, 3)). 
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In the CCC version of the SO(10) model , the mass matrices M u , Md, M e , Md, Ml, and 
Mr, for up quarks, down quarks, charged leptons, Dirac neutrinos, left-handed Majorana 
neutrinos, and right-handed Majorana neutrinos, respectively, are given, in the lowest tree 
level, by 

M u = SMQI) + A^m) + 1<0 3 + ')) + S^(</>%), 
M d = S^(^) + A^i-i^) + i(0 3 _')) " S^ 1 -^), 
rM e = S^(<P 1 _) + A( 12 °)(-(</> 3 _> - (0 3 _')) + S< 126 )(0 5 _), (23) 

r'M D = SW(<P 1 + ) + A^\(<Pl) - (0 s /)) - S^(<Pl), 
sM l = S^(<P(10,3,1)), 
s'M R = S^(4>(10,l,3)), 

where (0±) are the vacuum expectation values of the Higgs fields of 0io, (<f>±) and (0|') 
of 0i2o, and (0±), (0(10,3,1)) and (0(10,1,3)) of 0i 26 . See Ref [^] for details and the 



notations. The matrices S^ 10 ) and S^ 126 ) are real symmetric matrices and A^ 120 hs a pure 
imaginary matrix. These matrices are the coupling-constant matrices which appear in the 
Yukawa coupling of fermion fields with Higgs field , which is given by 

2L Y = 4 10) TjfofrMLi + 4 120) W^fa^Lj + S\f G) Wd c 4>i2^L 3 + (L^R)+ H.c. (24) 

The ipL(R)i are the 16 irreducible representations of the left- and right- handed fermion fields 
in a given i'th generation. The property that S^ 10 -* and 5 ( - 126 ^ are symmetric and A^ 120 ^ is 
antisymmetric results from the decomposition 16 x 16 = 10s + 120,4 + 126s, whereas the 
property that S^ 10 -* and S^ 126 ) are real and A^ 120 ^ is pure imaginary is a consequence of their 
being Hermitian, which in turn comes from the requirement of the invariance of Ly under 



the discrete symmetry ipi <-> iJ) R [20|. In Eq. (p3|) , the factors r ~ (2 ~ 3), r' , s and s' 



, all roughly of order unity, are the renormarization-group-equation factors [fM| [T^] which 
arise from the differences in the renormalization of the lepton and quark masses due to the 
color quantum numbers of the quarks and so on. The overall factor comes from the loop 

10 



correction of gauge boson in the renormalization group equation. Exactly we should consider 
the evolution equation of Yukawa coupling and in this case mass matrices gets renormalized 
in somewhat different form. Therefore, this form is an approximation. In this point we will 
also discuss in the last section. 

We now make the following assumptions. 

(i) The contribution from 120 is assumed to be small compared with the contributions 
from 10 and 126, and hence it is neglected in M u and Mb- On the other hand, it is retained 
in Md and M e , for the main term S^ 10 ' {<j>\) is smaller by the factor a = {(j)\) / (</>+) , which is 
of order of (rrib/mt) [see Eq.fl27|)]. This is an assumption for simplicity in order to incorporate 
Eq.® 

(ii) All the vacuum expectation values of Higgs fields are assumed to be real so that all 
the fermion mass matrices are Hermitian. 

With these assumptions, Eqs. (|24f) becomes 



M U = S + eS', 

M d = aS + S' + A' = aM u + A' - (ae - 1)5", 
rM e = aS- 35' + 6 A' = aM u + 6 A' - (ae + 3)5", (25) 

r'M D = S- 3eS", 
sM L = (3S\ 
s'M R = 1 S'. 

where the matrices S, S' and A' and the real parameters a, j3, 7 and 5 are defined by 



A' = A( 12 °)(-(0i) + ^<0 3 _'}), (26) 

« = (0i>/<0+>, P = (0(10,3, l))/(-i(0 5 _)), 

7 = (0(10, 1, 3)>/(-i(0 5 _», 6 = ((0 3 _) + m)/(m - \m). 



Note that solving diagonal elements of Eq . fl25|) for a, one finds 
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3TrM d + rTrM e _ m^ 
3TrM u + r'TrM D ~ m t ' 

which is about 0.02. As mentioned already, this is why the A^ 120 ' and S^ 12& > terms are kept 
in Md and M e . The Eq.(^) is our SO(10)-motivated model for fermion mass matrices. 

B. The case of 50(10) breaking down to SU(5) 

In this case, the fermion masses are also generated when the Higgs fields of 10, and 120, 
and 126 dimensional SO(10) representation (denoted by 0i O , 0i2O) an d 0i26 5 respectively) 
develop nonvanishing expectation values. Their decomposition under SU(5) are given by 

10 = 5 + 5, 
120 = 5 + 5 + 10 + 10 + 45 + 45, (28) 

126 = 1 + 5 + 10 + T5 + 45 + 50. 

The yukawa couplings in Ly gives the following fermion masses when the neutral components 
in a 5 and 5 Higgs multiplets in (fi w , 5, 5, 45, and 45 in 0i2o, and 1, 5, 15, and 45 in 0i 2 6 
of SU(5) C S*O(10) develop nonvanishing expectation values. 



S-j (ipLi) c 4>io*pLj -> S-j {{<f>io(5))(uR~iUL,3 + VrjVlj) + (0io (5)) (d R: id L ,j + eRje L ,j)}, 
A] {i>Li) c 4>i2oi)Lj -> A\f 0) {((j)i2o(5)}(d R ,id Ld + e^e LJ ) 

+((j>i2o(45)){d Rt idL,j - 3e^e L j) + (0i 2 o(5))z/^z/ Lj - 

Slj(i>Li) C 4>12(>1pLj, -> S , | J 126) {(0i26(5))(MR7ML,j - ?>T%iVL,j) 

+ (0126(45)) (dfl,idij - 3e^Ie Lj ) + (0126(1))^h,j*4,j 
+ (0 126 (15))^z/ L)i }, 

(29) 

Therefore, the mass matrices M u , M^, M e , Mo, M^, and Mr, for up quarks, down quarks, 
charged leptons, Dirac neutrinos, left-handed Majorana neutrinos, and right-handed Majo- 
rana neutrinos, respectively, are given by 
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M u = ^ 1O )(0 1O (5)) + A( 12 °)(0 12O (45)) + S( 126 )(0 126 (5)>, 
M d = S^(4> 10 (5)) + A( 12 °)((0 12O (5)> + (0 12O (45))) + S( 126 >(0 126 (45)>, 
rM e = S( 1O )(0 1O (5)) + A( 12 °)((0 12O (5)> - 3(0 12O (45)) - 35( 126 )(0 126 (45)), (30) 

r'M D = ^ 1O )(0 1O (5)) + A( 12 °)(0 12O (5)> -3^ 126 )(0 126 (5)>, 
S M L = 5 (126) (0 126 (15)), 
s'M R = S^(<p 126 (l)), 

These mass matrices reduce to the same form as Eq.fl2"5|) by assuming again that the con- 
tributions from 120 Higgs in M u and Mp are negligible and by defining the matrices S, S', 
and A' and the real parameters a, (3, 7, and 5 , instead of Eg. (|26|) , as 

S = S^(<p w (5)), S' = S( 126 )(0 126 (45)), 

A' = A( 12 °)((0 12O (5)) + (0 12O (45))), 

a = (<M5))/(0 1O (5)>, (3 = (0 126 (15)>/(0 126 (45)>, (31) 

7 =(0 126 (1))/(0 126 (45)), 
5 = ((0i2o(5)) - 3(0i 2O (45))/((0 12O (5)) + (0i 2O (45))), 

Thus, Eq.([25|) is our SO(10)-motivated model for fermion mass matrices both for the two 
SO (10) breaking patterns (i) and (ii) in Eq.([l9|). 



IV. FOUR TEXTURE ZERO MODEL IN SO(10) 



The SO(10) model Eq . (p25|) is now combined with the four texture zero ansatzae for M u , 
M d and M e which are given by Eq.(|6]). 

First it follows from Eq. (^5]) that S, S' and A' are represented in terms of the symmetric 
(antisymmetric) parts, M sym {M antisym ), of M u , M d and M e 



L ei 



(1 - ae)S = {M u ) sym - e(M d ) sym , 

S' = -{(M d ) sym - r(M e ) sym }, (32) 

I\ \^J-d)antisym- 
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We also find the constraints 

(1 - ae)r(M e ) sym = Aa{M u ) sym - (3 + ae)(M d ) sym , 

"\l'-ld)antisym '\^^e)antisym- \ / 

Using the four texture zero ansatzae for M u , Md and M e given by Eq.(BI), the respective 



elements of Eq . fl33[) become 



(1 — ae)rA e cos /3i 2 = 4aA u — (3 + ae)A d cos a± 2 , 

(1 — ae)rB e = 4aB u — (3 + ae)Bd, 
(1 — ae)rC e cos /3 2 3 = 4aC u — (3 + ae)Cd cos 0:23, 

(1 - ae)rD e = AaD u - (3 + ae)D d , (34) 

6 A d sin «i 2 = rA e sin /3 12 , 
<5C d sin «23 = rC e sin (3 2 3- 

In Eq.(|34|) there are six equations and eight unknown parameters, namely a, e, 5, «i 2 , a 2 3, 
Al2) /?23 and r provided that v4 u , B u ,...,D e are given. In the following, we treat cosai2 and 
r as free parameters so that all the other parameters are functions of them. Here we treat r 
as a free parameter too, although we know r ~ (2 ~ 3). Let us present the following useful 
expressions which are derived from Eq(^ 



(B e D d -D e B d \(A u \ (B e D u -D e B u \(A d > 

C ° SA2 = \B u Dd-D u B d ) W ~ \B u D d -D u B d ) U ) ""^ 
/ 5 e D d -D e B d \ (C u \ ( B e D u - D e B u \ ( C d \ 

COS/523 = \B u D d -D u B d ) [C-J ~ \B u D d -D u B d ) [c e ) C ° Sa23 ' 
sin/3 23 M e C d \ sin/5i2 



sin 0:23 V v4rfC e / sin an 



a = yv^ d -A^ d , _ ( __ — __ — _ () = r 



(35) 



\B u D d -D u B d ) 6 5= ( A e\ S m &2 

( bd db \ +1 > ' r( fg"-g^ ) ' " UJsina 12 - 

Now we discuss the MNS lepton mixing matrix and neutrino masses. The light Majorana 



neutrino mass matrix M v is given by Eq.(^), where the Dirac neutrino, left- and right- handed 
Majorana neutrino mass matrices M^, M^, and Mr are expressed in terms of the entries 
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of the quarks and charged lepton mass matrices due to the SO(10) constraints and their 
expressions are given, from Eqs.(|25|) and fl32|), by 

/ A D N 



r'M D = S- 3eS' 



A-d B D C D 
C D D D 



sM L = ps' = P 






A s , 





A s > 


B s > 


Cs> 





C31 


D s . j 



V 

/ A s > ^ 



s'M R = 7 S" = 7 



Ag/ B s > C s > 
C 5 ' As> 

where the elements A D , B D , Cd, Cd, A s >, B S ', Cs>, and Cs> are obtained as 



V 



/ 



(36) 



A D = A u - e(A d cosa 12 - rA e cos f3 12 ), 
B D = B u -e{B d -rB e ), 
C D = C u - e(C d cos «23 - rC e cos #23), 
D D = D u -e(D d -rD e ), 



(37) 



and 



B s > 
C s > 
D s > 



-(A d cosa 12 -rA e cos(3 12 ), 

\(B d -rB e ), 

- (C d cos a 23 - rC e cos (3 23 ) , 

\{D d -rD e ). 



(38) 



In the following analysis, we assume that the contribution of Mi to M v in Eq. (|3|) is 
much smaller than that of the second term so that we have M v = Mi — M^M^Md ~ 
—Mj)M^M D . Then, all the components A v , B v , C v , and D u in Eq.(|J) are determined, 
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from Eqs.([36|), fl37|), and (p8|), as functions of cos«i2 and r except for the common overall 
factor s'/(r /2 7) as. 



^ A v ^ 



M„ 



o a a 



- S '/(r' 2 7 ) 



7 

A D 
^4_d -Bd Cd 
C D D D 



A S / 

As" -^5' Cfi" 

C s , D s > 



^'0 A D 0^ 

Ad B d C d 
C D D D 



(39) 



A v 
B u 



' A 2 ' 



r /2 7 \ As'/ 
s' fA^Bn (C D A D C S > 



J-f2ry 



+A 



As' \Ds' As' Dgi J 

B D Cd Cs' A d ( — Cs> + B s > Dt 



D 



A S i A S i D s > As> D S ' 






s' \A d Cd 



rpt2ry 



A S ' 
i 

D 



D 



Cd Ad Cs' 



d 



D s > A S i D s > 



D 2 



(40) 



r /2 7 y Ds> J 

Therefore, the neutrino mass ratios mx/m^ and 1712/1713 and hence MNS lepton mixing 
matrix elements are also determined as functions of cos«i2 and r . The common overall 
factor s'/(r f2/ y) is determined by the Am 2 data from neutrino oscillation experiments. The 
light Majorana neutrino masses are obtained by diagonalizing M v as 



1 A v ^ 



O 



T 



Ay Dy Uj> 

a d v 



mi S 



O v 



J 



\ 



m 2 
m 3 



(41) 



For the case in which B U ,C U , D v ^> A v is satisfied, the neutrino masses are approximately 
expressed in terms of A u , B u , C u , and D v as 



m\ 



D v Ai 



B v D h 



c 



2' 
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m 2 ~ -{B v + D U - J{B V + D v )* - A(B U D U - Cl)} 



+ 



Dl + Cl - \D V {B V + D„ 



B v + D v y-A{B V D V -Cl)} 



{B V D V - Cl)j{B v + D v y-±{B v D v -Cl) 
1 



A 



(42) 



m 3 



{B v + D U + J(B U + D v f - A(B U D V - Cl)}. 



Dl + Cl - \D V {B U + D U + J{B P + D v y - A{B V D V - Cl)} 



Al. 



{B V D U - Cl)^{B u + D v f - A{B U D U - CD 
The orthogonal matrices O v which diagonalizes M v are expressed in terms of the diagonalized 



masses mi, m 2 , and 7713 and the matrix components A u , B v , C v , and Z)^ as [JT4 

/ 

— (( A.i-->, 

m 2 

(0,) 22 



0,. 



^(a) 2 i 
(a) 2 i 



ip v ] 



22 



A„ 



fe(^) 



\ 



in :i 
{O v ) W 



2:! 



with 



\ mi-D„ (°«')21 m2 _D„ ( ^)22 ma _",[,„ (0„) 



(a)L 



2.--! 



(43) 



(a)L 



(O 



W23 





1 






(— ) 2 


+ 1 + ( Cv 
1 


-) 2 


( — ) 2 


+ 1 + f C " 

1 


-) 2 



(44) 



'Al)2 1 J 1 c cv - 
^m.3' \mz—D U ' 



It should be remarked that the light neutrino mass matrix M v itself is out of type I via the 
seesaw mechanism and that the MNS lepton mixing matrix is obtained from Eqs. fl43|) and 
(|17|). Since O e is almost diagonal, the magnitudes of off-diagonal elements are predominated 
by Eq.(|44]). Thus the seesaw mechanism changes the form of lepton mixing matrix from that 
of CKM matrix given by Eg. (|i~6|) . 

Now, by changing the values of the free parameters in our model, we proceed to find the 
solutions which are consistent with the recent following findings that (i) the atmospheric 
neutrino oscillation experiment indicates the v [l -v r large mixing (0.28 < | C/23 1 ~ 0-72 p6[ ) 



with Am 23 



rrin 



rrin 



;i.5 ~ 6) x 10" 



experiments imply the MSW small mixing angle solution 



3.5 x 10 3 eV 2 , and (ii) the solar neutrino 

2 



271] with Aml 2 



m^ — m\ — (4 ~ 



10) x 10 6 eV 2 and sin 2 2#i 2 = (2 ~ 10) x 10 3 , or the large mixing angle solution [27]] with 
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Am 2 2 = m 2 — m\ ~ (8 ~ 30) x lCT 6 eV 2 and sin 2 26*i 2 = (0.5 ~ 1). In the following analysis, 
we transform A e , B e , C e , and D e in Eq.@ into — A e , —B e , —C e , and — D e , respectively by 
rephasing of the right-handed charged lepton fields. 

First assuming that the mass matrices M u , M^ and M e are all of type I, we calculate 
numerically the MNS lepton mixing matrix U using the central values for the running quarks 
and charged leptons masses at n — m z 



m u {m z ) = 2.33lg;^MeV, m c {m z ) = 677j|$MeV, m t {m z ) = 181 ± 13GeV, 
m d (m z ) = 4.69±{j;«JMeV, m s (m z ) = 93.4i^;gMeV, m b (m z ) = 3.00 ± O.llGeV, (45) 
m e {m z ) = 0.487MeV, m^(m z ) = 103MeV, m T {m z ) = 1.75GeV. 

Since the recent atmospheric neutrino oscillation data indicates large value of (2,3) element 
of U, (0.28 < \U 2 s\ < 0.72 [f26|j ), we obtain the allowed region of the parameters space, cosai 2 
vs r space which reproduces a large 1 1/23 1 - The result is given in Fig. 1. In this allowed 
parameter region, r ~ 3 is automatically satisfied without any fine tuning. However, we 
have a serious problem that in this allowed parameter space we cannot accommodate the 
overall factor s' /(r' 2, y) simultaneously to the data Am 2 2 = ^2 ~~ m i ~ 10~ 4 eV 2 (Here we 
have adopted a rather conservative value. We accept more restrictive ones later.) from 
solar neutrino oscillation experiments and the data Am| 3 = rnj — m 2 , ~ 3.5 x 10 _3 eV 2 from 
atmospheric neutrino oscillation experiments. Taking deviations from the central values 
p8| for quarks and charged leptons masses does not resolve this problem. This difficulty is 
resolved by abandoning the above assumption that the mass matrices M u , Md and M e are 
all of type I. Let us assume that M e deviates from type I although M u and M d are of type 
I. Then, we can accommodate the overall factor s'/(r' 2/ ~f) simultaneously to both Am 2 data 
from solar and atmospheric neutrino oscillation experiments. 

Next we discuss this new scenario and show that there are solutions consistent with the 
data. First we represent the deviation from type I as B e = m M (l + £). In this case, the 
entries of the mass matrix M e for charged leptons in Eq.@ are given, in the unit of eV, as 



/ m P m,,m T R 

A e = J e -^^ ~ 7.1 x 10 6 , 



B e = (l + £)m„ ~ 1.02 x 10 8 (1 + 

a 



m T — £ m. 



, (m e + f m M ) (m r - (1 + £) m M - m e ) I ""\ s "" M ~ V 8 - 1 x 1Ql4 + L7 x 1Ql7 ^ 

£) e = -m e -£ m^ + m T ~ 7.1 x 10 9 , (46) 

and the expressions for e, a and the elements of the matrices S and 5" are given by 

(— r(l + £)m^ + 3m s )m t — m c (3m b — rr7i r ) 
r((l + i)m b m^ - m s m T ) 
r((l + C)m b m^ - (1 + £) m d"V _ m s m T ) 



a 



m b m c + (—rrrifj, — r^m^ — m s )m t — m c {m d — rm T ) ' 
^ — ((rm^ + m s )(cosfti 2v /m^mI(-3mbm c - r(l + Qm^rnt) + 

cosfti2A/ m <^s(3m s m t + rm c m T ) — r(—m.bm,j, + m s m T ) y/m c m u )) / 

{Arm s {m b m^ — m s m T )), 



S- 



22 ^ 



{—rmp - m s )(-r(l + £)m^ + 3m s )m t 



4r(mbm M — m s m T ) 
S23 — ( - cos a 2 3^m b m d (-rm ll - m s )m t 

(—3m b m c — r(l + £)ra M mt + 3m s m t + rm c m T ) — 

r(-rm^ - m s )m t ((l + Cjm b m^ - m s m T )^/m t m u )/ 

{r 2 m T {rn b m c m iL + m s {^m^m t — m c m T )) — 

r{3m b 2 m c m^ + Am s 2 m t m T — m b m s ((A + Qm^mt + 3m c m T ))), 

{—rrrifj, — m s )mt(3m b — rm T ) 
33 '" 



5»»~ v ' M 



4r(m( ) m /i — m s m T ) 



(cosfti 2v /m d m s (rm M + n^mj - r{m b m ii - m s m T )y/m c m u ) 

'-'12 ~ ^ 

4m s m( 



1 

4 l 

r{m b m^ — m s m T )m v 



S'22 — -( rm n + m 



^ - _ COSft 23V /m 6 m d + / . 

V (-rm M -m s )y/m t m u 

( (m b m^ - m s m T )((-r(l + £)m^ + 3m s )m t - m c (3m b - rm T )) 

V (-rm M -m s )m t ((l +£)m fe m At -m s m T ) 

£33 - -(m b -m d + rm T ). (47) 



The point r ~ 3 is rather singular in the following sense. As is seen from Eq.fl3"5D, e becomes 
small at r ~ 3, hence we can not neglect £ in this region. That is, e is sensitive to small £. 
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For instance, substituting Eq.(f47|) into Eqs.(p5"D we obtain Mo and Mr and therefore M v 
through Eg. (|39"D . The behaviors of the elements in the neutrino mass matrix M v are depicted 
in Fig.2. The large v^-v T neutrino mixing appears under the condition that B v ~ D v which 
is realized at £ ~ 0.02. By changing £ freely with fixed r(= 3), we can well reproduce the 
experimental data as shown in Figs. 3-5, in which the constraints from | C/23 1 5 A m i2/^ m 23! 
and the both are satisfied, respectively. It is seen from Fig. 5 that by deviating B e a little bit 
from type I (£ ~ 0.01), we can well reproduce the experimental data for the solar neutrino 
oscillation and atmospheric neutrino deficit. If we relax the condition r = 3 and change r 
freely as well as cos a.12 and £ around the values of the above solutions, we have the larger 
allowed region as is shown in Fig. 6. In the above allowed regions shown in Figs. 1-5, we have 
used only the conservative condition for Am^ from the solar neutrino experiments, that is, 
we have not used the constraints of the mixing angle from the solar neutrino oscillation 
experiments. When we take them into account in addition to the constraints from Am^, we 
obtain more restrictive allowed region than that of Fig. 6. Under the condition of the small 
mixing angle solution for solar neutrino experiments, the larger region of | cosa^l in Fig. 6 
is eliminated and we have the allowed region as is shown in Fig. 7. On the other hand, under 
the condition of the large mixing angle solution, the smaller region of | cosa^l is eliminated 
and the allowed region is given in Fig. 8. It should be noted that as seen in Figs. 6-8 our 
model not only satisfies the experimental observations in the lepton sector but also provides 
the restriction on the CP violation phase, cos«i2, from the neutrino oscillation experiments. 
Of course, we can also restrict the other CP violation phases, cosa 2 3, cos/?i 2 and cos/3 2 3, 
which are respectively depicted in Fig. 9-11. Also it follows from Eq. ([39]) that the neutrino 
mass ratios \m\jm2\ and \m2/m^\ become sensitive functions of £, as are shown in Fig.12 
taking typical values of r and cos«i2. The common overall factor s'/(r' 2/ ~f) in Eq.(|39"D is 
determined to be of order 10~ 13 by the Am 2 data from the solar and atmospheric neutrino 
oscillation experiments. 

Finally we discuss the entries of the CKM quark mixing matrix which are given by 
Eqs.fllB]) with (J35l). In our model, all the elements of the CKM mixing matrix are also 
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functions of two free parameters cos ol\i and £. The parameters determined so far from the 
lepton sector do not give rise any inconsistency with the data in quark sector. 

V. SUMMARY 

In this paper we have presented and discussed a model of texture four zero quark-lepton 
mass matrices in the context of SO (10). The consistent fitting of the free parameters to 
the data for neutrino oscillation experiments has forced us to use the charged lepton mass 
matrix which slightly deviates from purely type I form (£ ~ 0.01). Using this deviated type 
of mass matrix for the charged leptons and the mass matrices for quarks of type I, we have 
been able to reproduce four entries in the CKM quark mixing matrix and to predict six 
entries in the MNS lepton mixing matrix and three Majorana neutrino masses which are 
consistent with the experimental data. The model has also given the restrictions on the CP 
violating phases which came from the neutrino oscillation experiments. Remarkably enough 
the parameter r fixed from data fitting is coincident with the value r ~ (2 ~ 3) obtained 



from the renormalization equation [Q . So it is attractive to expect that the above deviation 
(£ ~ 0.01) from type I form can be obtained by taking the evolution equation of Yukawa 
coupling fully. (In this paper we have considered the loop correction of gauge boson in the 
evolution equation.) Though the detail calculations will be developed in the forthcoming 
paper, we will roughly outline our idea. That is, (charged lepton) mass matrix is exactly of 
type I at some scale. However, they change their form due to the evolutionary equation of 
the Yukawa coupling Y a until the corresponding Higgs field acquires the vacuum expectation 
value 



where T*, G*, H' are the vertex corrections due to the fermion, the gauge boson and the 
Higgs boson, respectively. After that, each mass furthermore changes its value according 
to the mass renormalization equation. The evolution equation of Yukawa coupling is very 
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sensitive to the Higgs potentials and the initial conditions. One such sensitivity has been 
found in the behavior of £. The detail will be given in the forthcoming paper. 
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FIGURES 
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o > ■ ■ ■ ■ ■ ■ ■ 1 

-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 



COS 0|2 

FIG. 1. The case where M u , M^, and M e are all of purely type I is analyzed. The experimental 
constraint on | C/23 1 ( 0.28 < 1 1/^23 1 2 < 0.72 ) gives the allowed region (shaded area) in the cos«i2-r 
plane. Here the r is treated as a free parameter. 
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FIG. 2. The slight deviation from type I (£ 7^ 0) makes physical parameters change drastically . 
The dotted lines (solid lines) show the £ dependence for cos «23 > (cos «23 < 0) in each diagram. 
All lines terminate at the points from where | cos 023 1 > 1 or | cos /?23| > 1 as will be seen from Fig. 9 
and Fig. 11. (a) The diagram of the elements in the neutrino mass matrix M v versus £. Except for 
\C u (r' 2 j/s')\ the dotted lines are overlapped with the corresponding solid lines, (b) The diagram 
of the neutrino mass eigen values versus £. (c) The MNS mixing matrices versus £. 
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FIG. 3. The experimental constraint on IL^I ( 0.28 < IC/23I 2 < 0.72 ) gives the allowed region 
(dotted area) in the cosai2-£ plane. 
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FIG. 4. The allowed region in the cos ai2 - £ plane from the experimental constraints 
Am? 2 /Am| 3 < (1 X 10~ 4 )/(3.5 x 10" 3 ) = 2.9 x 10~ 2 . 
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FIG. 5. The allowed region in the cos«i2 - £ plane from the experimental constraints 
0.28 < \U 2 3\ 2 < 0.72 and Amf 2 /Am| 3 < 2.9 x 10~ 2 . 
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FIG. 6. The r is treated as a free parameter, (a) shows the allowed region in the cosai2 - r 
- £ space from the experimental constraints 0.28 < 1 1/^23 1 2 < 0.72 and Am 2 2 /Am| 3 < 2.9 x 10 -2 . 
(b), (c) and (d) show the projected allowed regions in the cosai2-£, £-r, and cosai2-r planes, 
respectively. 



(b) 

0.2 




uj 0.1 



(d) 



0.4 -0.2 0.2 0.4 

COS a 12 



H 3 



2.5 



0.2 -0.4 -0.2 0.2 0.4 

COS a 12 



FIG. 7. The allowed region in the cosai2-r-£ space from the experimental constraints 
0.28 < 1 1/23 1 2 < 0.72, the small mixing angle solution of the solar neutrino experiments 
(sin 2 2^12 = (2 ~ 10) x 10~ 3 ), and the up-to-date value of mass difference Am^/Am^ = 
((4 ~ 10) x 10- 6 )/((1.5 ~ 6) x 10~ 3 ) = (0.67 ~ 6.7) x 10~ 3 . 
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FIG. 8. The allowed region in the cosai2-r-£ space from the experimental con- 

straints 0.28 < 1 1/^23 1 2 < 0.72, the large mixing angle solution of the solar neutrino exper- 
iments sin 2 20i2 = (0.5 ~ 1), and the up-to-date value of mass difference Am^/Am 2 ^ = 
((8 ~ 30) x 10- 6 )/((1.5 ~ 6) x 10~ 3 ) = (0.13 ~ 2.0) x lO" 3 . 
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FIG. 9. The allowed region in the cos«23 - r - £ space from the same constraints as in Fig. 6. 
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FIG. 10. The allowed region in the cos/3i2 - r - £ space from the same constraints as in Fig. 6. 
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FIG. 11. The allowed region in the cos/323 - r - £ space from the same constraints as in Fig. 6. 
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FIG. 12. The dependence of log 10 \m\/m2\ (a) and log 10 \m2/m^\ (b) on £ for cosai2 = 0.2 
and r = 3. The dotted line (solid line) shows the £ dependence for cos 023 > (cos 023 < 0) in 
each diagram. Both lines terminate at the points from where | cos 023 1 > 1 or | cos /323 1 > 1 as are 
seen from Fig. 9 and Fig.ll. The singular behaviors of (a) and (b) come from those of m\ and 777,2 
(see Fig.2). 
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